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ABSTRACT 

Current precision in radial velocity (RV) measurements of binary stars reaches ~2 ms"'. This 
level of precision means that RV models have to take into account additional non-Keplerian 
effects such as tidal and rotational distortion of the components of a binary star, relativistic 
effects and orbital precession. Such an approach is necessary when one wants to search for 
planets or precisely measure fundamental parameters of stars with a very high accuracy using 
precision RVs of binary stars. 

We generate synthetic binaries using Yonsei-Yale stellar models. For typical represen- 
tatives we investigate the impact of various orbital orientations and different non-Keplerian 
effects on the RV curves. To this end we simulate RV observations with an added white noise 
of different scale. Subsequently we try to reconstruct the input orbital parameters and their 
errors by fitting a model using a standard least-squares method. In particular we investigate 
the connection between the tidal distortion of the shape of the stars and the best-fit orbital 
eccentricity, the possibility of deriving orbital inclination of a non-eclipsing binary star by 
exploiting relativistic effects and the circumstances in which the orbital precession can be 
detected. 

We confirm that the method proposed by IZucker & Alexander! (1200 7h to obtain orbital 
inclination with use of the relativistic effect does work in favourable cases and that it can 
be used even for orbital configurations far from an edge-on orientation. We show that the 
RV variations imposed by tidally distorted stars can mimic non-zero eccentricity in some 
binaries. The scale of such an effect depends on the RV accuracy. Finally, we demonstrate that 
the apsidal precession can be easily detected with precision RVs. In particular we can detect 
orbital precession of 10 '* rad yr"', 10"^ rad yr ' for precision of RVs of 1 ms"' and 10 ms"' 
respectively. 

Key words: binaries: spectroscopic - stars: fundamental parameters - relativity - methods: 
analytical - methods: numerical - techniques: radial velocities. 



1 INTRODUCTION 

In 1842, Christian Andreas Doppler published his work "On the 
colo ured light of the binary stars and some other stars of the heav- 
ens" jPoppleil 18421) . Although he did not get it quite right then, his 
work has started the career of the well known Doppler effect. Over 
a century has passed and the best radial velocity (RV) precision for 
the primaries of double-lined spectroscopic binaries remained al- 
most unchanged at about 100 ms"' until recently. In 2003, Zucker 
et al. presented a precise radial velocity measurement of a binary 
component with precision reaching 10 ms"' . With the use of TOD- 
COR, high quality echelle spectra from CORALIE and faint second 
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component not interfering with the lines of primary component the 
new capability was reached. A year later Skuljan et al (2004) have 
reached a precision of 14 ms"' in a similar configuration. Another 
important step was development of spectra disentangling code by 
M. Konacki utilizing the iodine cell technique and making the reg- 
ular, high precision RV measurement of many binaries, with two 
similar components, possible. It is now within the reach of many 
instruments to achieve the precision of up to a few meters per sec- 
ond for the compone nts of double-lined spectroscopic binary stars 
jKonacki et allUOK]!) . 

The main goal of this paper is to investigate the selected fea- 
tures of double-lined binary stars, their potential and problems in 
the context of precise RV measurement s. We based our sele ction 
of effects on the recent publications bv lKonacki etalj | |201(]|) and 
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iMazed ( l2008h and observational practices. In the new precision 
regime the tidal effects become important. Stars cannot be treated 
as single points in space. The presence of another body in an orbit 
around a star generates tidal effects influencing the shape of a star. 
Because of it, the measured radial vel ocity is different from that de- 
rived f rom simple Keplerian motion jWilson & Sofia|[l976l : lKoD^ 
Il980ah . This effect, when neglected, may be responsible for the 
measurement of a non-zero eccentricity in an otherwise circular or- 
bit. Given sufficient precision of the RV curve, one can rule out this 
situation to a certain extent. Moreover the relativistic effects may 
be a source of additional informat ion about a system. T herefore, 
we explore the relations derived bv lZucker & Alexander] (12007') re- 
lating the relativistic effects and the orbital inclination, i. With suf- 
ficient precision of the RV measurements, one can derive ; of the 
system even for a non-eclipsing systems and using only RVs. The 
last effect we examine is the apsidal precession. We determine the 
required RV precision to measure the rate of apsidal precession. To 
achieve these goals we conduct thorough numerical simulations us- 
ing currently available codes for modelling binary stars. Our anal- 
ysis is aimed to support precise, current and future, RV surveys of 
binary stars by investigating the challenges and possibilities. 

In section 2, we describe models as well as a synthetic popu- 
lation of binary stars which are used in the tests. This section also 
explains an algorithm employed to obtain a final set of parame- 
ters based on artificial data. In Section 3, we analyse the impact of 
tidal effects on the R V curve and ecc entricity. In Section 4, we test 
the idea proposed bv lZucker & Alexan der (2007) of using pure RV 
data as an independent source of the orbital inclination. We test its 
feasibility and usability. In Section 5 we check the limits of apsidal 
motion detection for different RV precision. The conclusions are 
provided in Section 6. 



2 MODELS 

The tidal and rotational e ffects were s tudied and implemented as a 
computer algorithm by iWilsonI ( 19931) . The analytical approxima- 
tion was provided bv lKopall (1980a). We decided to use WD code 
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based on the models of Wils on & DevirmevI jl97 1! ) and Wilson 

( Il979h as the source of synthetic RV data. The WD code includes 
all the effects we are interested in and much more. It is moreover a 
widely accepted tool for modelling eclipsing binary stars. 

We prepared the artificial population of binary stars in differ- 
ent orbital configurations, containing from 10000 to 400000 ob- 
jects. The resolution of such a grid and the choice of parameters 
depended on an analysed problem and a computing time neces- 
sary to pass the grid. Every problem had been computed several 
times to confirm the results. The input para meters are based on the 
widely used Y onsei-Yale {Y^) isochrones (Demarque et al."2004'; 
Yi et al.ll200lh. We tested other ones li ke PADOVA (iGirardi et al.i 
200Cl;lMarigo et alj2008l) or GENEVA dLeieune & Schaerej200l'l) 
but the resulting differences had minor impact on the results and 
none on the conclusions. Chemical abundances in our set are de- 
scribed by X = 0.71, Y = 0.27, Z = 0.02, this composition is 
similar to the Sun. Certain values for the test stars were obtained 



by a linear interpolation of neighbouring values from the ta- 
ble. Additional information necessary to describe the components 
of a binary , like the gravi t ationa l brigh t ening coefficien ts, were 
taken from lAlencar & Va3 (ll997L Il999[) : Iciaretl ( l2000alibl) . Stars 
and their orbital parameters are described in Table [r|and|2]respec- 
tively. Limb darkening was compu ted with cosine law with WD 
code dWilson & Van HammdfioOTl) . Stellar atmosphere formula- 



tion was used for local emission on binary star components, in- 
stead of a black body approximation. Other internal parameters of 
WD code used in our simulations are summarised in Table |3] We 
assume constant F, periastron-synchronised ratio of the axial rota- 
tion rate to the mean orbital rate. The sample aims to explore the 
set of close binaries with masses between 0.5 Mq and 1.5 Mq. In 
this range the examined effects are important for precise RV mea- 
surements. Symbols used in the tables are the orbital elements, e 
for the eccentricity, uj for the argument of pericenter. The parame- 
ter ctrv describes the white noise added to every measurement and 
i is given in degrees. In the case of an eclipsing configuration, we 
chose orbits close to an edge-on orientation, i ranges from 87° to 
93°. This range assures that every object has an eclipse. For a given 
range of orbital inclinations (eclipsing or not), the configurations 
were drawn from a uniform distribution of the orbital inclination. 
This distribution differs from tlie isotropic distribution of pos- 
sible orbits in space. The results are not a representative statis- 
tical sample of real binaries, but rather demonstrate the kind 
of effects we expect to see in the sampled range of parameters. 

In our simulations we proceeded as follows. We generated a 
synthetic population of binaries and for each binary we computed 
a synthetic RV set with white noise added to mimic a measurement 
error. For every binary, we computed thirty RVs for the primary and 
secondary components randomly spread over a given time span of a 
data set. Such an RV set was used to fi t a model. We used the MIN- 
PACK library (More et al.ll980l[l984H'l to carry out a least-squares 
fitting. As input parameters we used the original ones modified by 
at least half of a percent. 

In the models we assume that the stars are separated and none 
of the components fills its Roche lobe. This is accomplished by 
testing the WD code output and fulfilling the inequality: 



2(/?i +^^2) <a(i-e) 



(1) 



where Ri, Ri are the radii of the first and second star respectively. 
This assures that we deal only with detached systems. We analyse 
RV curves for both stars at the same time. Moments of eclipses are 
avoided. It is a common approach made by observers which makes 
the RV modelling much easier. However, it is worth noticing that 
eclipses may be a reliable source of additional information espe- 
cially in the case of rapidly rotating non-circular binaries where 
the non-Keplerian effects have a high amplitude. This is known a s 
the Rossiter-McLaughlin effect ( lRossiteJl924l : lMcLaughlirJl924t) . 
Unfortunately WD code does not allow to model a misalignment 
between the star's spin and the orbital angular momentum of the 
binary, and covers only the basic case of aligned rotation. In our 
simulations we used as an input only the aligned cases. Therefore 
our results are not influenced by this limitation of the WD code. 
However during the normal analysis of observational data it should 
be included. Proper modelling of spectral lines and radial veloc- 
ity should take this effect into account, especially in the case of 
eclipses, where it is essential. Omitting it may cause the observed 
precession rate to be four times slower than the theoretical rate 
( Albrecht et al. 2009). The misalignment will change the shape of 
the RV curve but it's highest amplitude will remain the same. Mis- 
aligned configurations produce generally lower amplitudes. 

A nalytica l approach to the RV curves was described by 

[Kopal' l'l980a"P). More recent publications regarding the subject 
by Zucker & Alexander (2007) and Mazeh (2008) prove the signif- 
icance of modelling close interactions in binaries. Development of 
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Table 2. Orbital and other parameters used in the simulations. 



Varied Parameter Values Range Unit 



0.0, 0.1,.. ,0.9 - 
0.0 - 2jt 

1, 10, 100 ms-' 
3, 5, 10 day 
87, 88 .. 93 
0, 3 .. 90 
87 



P 

i (eclipsing binaries) 
! (relativistic eft'ects) 
I (apsidal precession) 



Table 3. Internal parameters of Wilson-Devinney code used in the simula- 
tions. 



Parameter 


Value 


Grid size for star 1 and 2 


30 


Bolometric albedo of star 1 and 2 


0.5 


F 





high res olution spectrogr aphs hke HARPS dRupprecht at al.l2004|>. 
HIRE S l lVogt et alJirggj) or the incoming ESPRESSO (Pepe et al. 

and methods like iodine cell calibration (■Konacki ,2005') en- 
sure the advancement of precision in RV measurements. Software 
and model evolution is also needed to fully utilise hardware im- 
provements. One exa mple of such an improvement is the spectral 
disentangling method ( iHadravall 19951 : iKonacki et allllOld^ . 



3 TIDAL EFFECT AND ECCENTRICITY 

The shape of stars in close binaries differs from standalone objects. 
Both components are subject to mutual tidal forces and also mutual 
irradiation. Deviation from the spherical shape and non-uniform 
surface brightness are a source of an RV curve distortion. An av- 
erage amplitude of the tidal effect in our sample is 174 ms"', but 
can be as big as 3800 ms"' or as small as 0.25 ms"'. We define 
tidal effect semi-amplitude as the maximal difference between the 
Keplerian RV and the RV curve of the same star with the tidally 
distorted surface. In our simulations this parameter is estimated on 
the base of randomly distributed measurements across the system's 
period. Examples of three cases are presented in Figure [T] panels 
(b-d) respectively. When this effect is ignored one can obtain a best- 
fit RV root mean square (RMS) of 1 ms"' or better, only in the case 
of well separated objects. 

Moreover, the tidal effect may mimic an eccentric orbit when 
a simple Keplerian model is fitted to the RV data. The tidal effect 
is often easily incorporated into the eccentricity. We show to what 
extent the tidal effect may produce correct fits with RMS of the 
best-fit (Ty„ lower than Sctrv ■ To this end, almost 12000 systems had 
been tested. The synthetic population consists of all age and mass 
permutations from Table [T] used with the entire range of varied 
orbital parameters from Table |2] Only the eccentricity e is kept at 
zero. 

From our simulations emerges a safe limit for a separation of 
stars: a(l-e) > 20(R[ +R2)- In this regime, we may safely ignore the 
tidal effects and assume them to contribute less than 0.1 ms"'. The 
tides' impact on the RVs starts to drop rapidly with the parameter i- 
lower than 0.1, where s is given by the equation: 
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Figure 2. Illustration of the relation between the semi-amplitude of the tidal 
effect and the « parameter. Almost 12000 simulations in this figure let us 
estimate the safe limits for avoiding the tidal eftect at a given RV precision. 



R, +R2 
a{l - e) 



(2) 



A detailed relation between the s parameter and the semi-amplitude 
of the tidal effect is illustrated in Figure |2] 

Figure [3] shows the results of simulations and three cases 
where the best-fit RMS, o"/,,, is within a reasonable limit of 2<rf,v 
but the fit still returns a non-zero eccentricity. These three cases. 
Figure [3] (b, d, f), illustrates the limitation of a simple Keplerian 
model, applied to the tidally distorted stars. We also used another 
marker to rule out most of the synthetic binaries with e > 0: 
e/cTe > 3. The error of the fitted eccentricity cr^ is very often high 
enough to still accommodate zero as a value within icr^ limit. How- 
ever, for a-Rv = 1, 10 and 100 ms"', 28%, 32% and 10% of the 
best-fit e, respectively do not meet this criterion. This leaves us with 
incorrect results that may be treated as a detection of an eccentric 
orbit. 

We propose to use three safe limits for ruling out eccentric 
false-positives. They are computed as three times the highest false- 
positive e from our sample. It is e = 0.0004, 0.003 and 0.01 for the 
three RV measurements precisions 1, 10 and 100 ms"'. The limits 
are illustrated in Figure |3](b, d, f) via horizontal green lines. One 
may linearly interpolate given limits, but we do not recommend ex- 
trapolation as the simulations were conducted for a specific set of 
binaries. An eccentricity lower than this limit and even complying 
with all the restrictions, needs special attention and more advanced 
model than a simple Keplerian one should be used. Further statisti- 
cal tests (e.g. ILucvI ( l2005b ) may be used to decide if a more sophis- 
ticated model is necessary. It is worth noting that with light curves 
of eclipsing binaries one may directly measure the eccentricity and 
argument of periastron. This procedure requires times of minima 
for both components. Such an eccentricity may be compared with 
the one derived from RVs. 



4 RELATIVISTIC EFFECTS AND ORBITAL 
INCLINATION 

In the basic Keplerian model, the equations for the RV of the pri- 
mary star Vzi and the secondary star RV Vzi, in the direction z, 
are: 



Vzi,2 = ^i.2(cos (tJi.2 -I- v) -I- ecos 0)12) + Vzi 



(3) 
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Table 1. Stellar parameters from Y- models for the components of synthetic binaries. Permutations used in our simulation assume the same age of components. 



Parameter Star 1 Star 2 Star 3 Star 4 Star 5 Star 6 Star 7 Star 8 Unit 



Age 


0.8 


0.8 


0.8 


3 


3 


3 


12 


12 


Gyr 


Mass 


0.5 


1 


1.5 


0.5 


1 


1.5 


0.5 


1 


Mo 


Logg 


4.78 


4.53 


4.22 


4.77 


4.49 


3.34 


4.75 


3.68 




Temperature 


3521 


5619 


6995 


3518 


4927 


5676 


3563 


4823 


K 


Gravitational brightening 


0.12 


0.396 


0.272 


0.12 


0.396 


0.4 


0.12 


0.38 





[ht] 






0.25 




0.2 




0.15 




0.1 




0.05 


to 




1 





> 






-0.05 




-0.1 




-0.15 




-0.2 




-0.25 



upy = 1 m/s 
f 1 m/s 
Orv =i 100 m/s 
const i = 3 m/s 
const 2 ^ 30 m/s 
const 3 =i 300 m/s 



10000 



RMS of the fit [m/s] 
(a) 



M,=0.5 Mgun 
M2=0.5 Ms„„ 



if 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
phase [-] 

(c) 




4000 
3000 
2000 
1000 


-1000 
-2000 
-3000 
-4000 



^2=1 Msun 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
phase [-] 

(d) 



Figure 1. Illustration of the tidal effect's impact on fitting just a simple Keplerian model to the RVs of a binary star. The panel (a) shows the full result of 
our simulation, a relation between the best-fit RMS of the final fit and the eccentricity. Synthetic data sets are based on WD code while a simple Kepleiian 
model is used for modelling. Vertical lines denote three constant RMS intervals 3, 30 and 300 ms"'. All the e results on the left side of each line may be easily 
taken as a genuine eccentricity. The size of each dot is determined by the RMS of the fit. The panels (b-d) show three typical examples of tidal effect in RV: 
average, minimal and maximal. As one may easily predict, more massive, tight and old binaries tend to have much bigger surface distortions, and because of 
this, higher semi-amplitudes of the tidal eft'ect. Gaps in the curves are caused by eclipses where no measurements were modelled. The system shown in the 
panel (c) has two shallow and short eclipses around phase 0.3 and 0.95, due to the orbital configuration and small size of the components. Panel (b) shows 
additional information, scaled down by factor of 1000 RV coming from Keplerian motion of components around common center of mass. 



where radial velocity semi-amplitudes Ki and are equal to: 



P{q + 1) Vl -e^ 

Ki = qKi (5) 

a is the semi-major axis of the component's orbit, q mass ratio of 
the binary, P stands for period and e is the system's eccentricity. In 
order to fully describe the system we need the argument of perias- 



tron a>, the radial velocity Vzo of the center of mass and the true 
anomaly v. 

The tidal effect is not the only subtle effect that has to be 
considered while modelling precision RVs of binaries. Relativis- 
tic effects come into play often at the level well above the 1 ms"'. 
At least three relativistic contributions have to be considered. The 
light-time (LT) effect within an orbit and coming from the finite 
speed of light: 
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(d) 



0.01 
0.001 
0.0001 
1e-005 
1e-006 
1e-007 
1e-008 
1e-009 
1e-010 




o o 



o 
o 



■ RMS = 1 m/s 
o RMS = 5 m/s 
o RMS = 1 m/s 
O RMS = 1 00 m/s 
O RMS = 1 000 m/s 



-3.2 



-3.4 



-3.6 



-3.8 



0.1 1 10 100 1000 10000 

RV semi-amplitude of tide effect [m/s] 

(e) 



0.01 
0.001 
0.0001 
1e-005 
1e-006 
1e-007 
1e-008 
1e-009 
1e-010 




Orv = 1 00 m/s 
const = 0.01 



50 60 70 80 100 200 
RMS of the fit [m/s] 

(f) 



Figure 3. A more detailed view of the simulations of the tidal effect. The panels (a, c, e) present all the results from the synthetic set for the three different 
values of white noise, ctrv = 1,10 and 100 ms"' respectively. Each panel shows the relation between the highest out of two components semi-amplitude of the 
tidal effect and the best-fit eccentricity. The corresponding figures in the right column are the results of the simulation with 2itrv limit applied. The size of the 
symbol in the first column is defined by the best-fit RMS given by log RIViS . The size of the symbol in the second column, panels (b, d, f), is constant. Colour 
denotes the error of the best-fit eccentricity <Tg as log a-g . Results clearly show that the eccentricity may be incorporated into the simple model and no hint of 
an inadequate model is visible. With the increasing semi-amplitude of the tidal eft'ect, increases the best-fit RMS, <Tfi,. One may notice that the semi-amplitude 
even a few times bigger than the white noise ctrv, returns a reasonable results with low RMS and low cr^. It may be hard to distinguish errors coming from the 
underestimated RV errors and a too simplistic model. It is worth noticing that as the best-fit RMS starts to grow the error of best-fit eccentricity decreases. 
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Kf2 2 

AVt7-i2 = — ^ sin (v + Wi 2)(1 + e cos v) (6) 
c 

This is a correct approximation to the order, where fl = v/c. The 
transverse Doppler effect: 



AVt 



KI2 /, l-e^ 
1 + ecos y — 



and the gravitational redshift: 

. , , ^l,2(2A'i.2 + A^2,l) 

AVg«i_2 = — ^ e 



(7) 



(8) 



The speed of light is denoted by c. On e may find m o re de- 
tai ls on those equation s in th e papers by iKonacki et al] ilOld) 
or IZucker & Alexander! ( l2007h . The principals of the relativis- 
tic effects on the motion o f a binary star are well described by 
iKopeikin & Ozemo\i ( Il999l) . We omit other relativistic effects of 
orders higher than jS^. 

We have used Equations (|6} - ((Sj to generate a set of RVs of 
over 400000 synthetic close binaries and explored the average scale 
and impact of the above effects on the RV curves. In our illustra- 
tion we chose two systems with different eccentricities. In the first 
case where only LT effect is visible e = and in the second case 
e = 0.2. In the first case TD and GR effects are constant but in the 
second case vary with time. Other parameters as well as images il- 
lustrating the contribution of each effect to the final RV are shown 
in Figure |4] Our main interest is however in a way to obtain the 
inclination of an orbit purel y from RV measurements, via equation 
jZucker & AlexandedlioOTh : 



3e 



where 



— arctan 



(9) 



(10) 



The new elements Kg i,Ks2,Kci and Kc2 are the modified elements 
Ki and K2 from the simple RV Equation [3] now accommodating 
relativistic effects. A detailed description and in terpretation of these 
new elements is in Zucker & AlexandeJ(l2007h . 

In these simulations we assumed the three following condi- 
tions: 

• The total RMS of fitted relativistic model is lower than three 
times the white noise applied to the RV data: cr^„ < Actrv, where 
A =3. 

• The error of inclination tr, is ten times lower than inclination 
;: ; > Scr,, where B = 10. 

• The difference between the true inclination /, and the best-fit 
inclination is lower than 3o"/: |j, - ;| < Ccr,, where C =3. 

The result of such a simulation is presented in Figure |5] It clearly 
shows that the proposed method works, but obviously only for ec- 
centric orbits. The precision of the fit can be as good as 0.01° for 
1 ms"' accuracy of the RV data and 0.1° for 10 ms"' accuracy. For 
100 ms"' precision there is no single system in the final pool. How- 
ever these values are possible only for very low inclinations, lower 
than 20 degrees. Such systems will have spectral lines blended, 
making the measurement of RV a very hard or even impossible 
task. Altogether it means that 40% and 25% of results of respective 
1 ms"' and 10 ms"' precisions meet the conditions. We are in a 
good situation because we know the true inclination. However re- 
moving the third condition does not change the whole picture and 
adds only a fraction of percent to the results. When we loosen our 



Table 4. Tfie impact of different tliresliolds on the number of acceptable 
solutions for the orbital inclination. The columns 5-7 represent a percent of 
systems from the synthetic set where all thi'ee conditions with respective A, 
B and C values were met. The typical number of systems per one ctrv value 
and per one set is 51000, where 10% of them have a circular orbit. The 
relativistic effect is much more sensitive to small inclinations than to large 
(edge-on orbits) inclinations. The effect strongly depends on the ability to 



Paiameter Set 


A 


B 


C 


1 ms"' 


10 ms"' 


100 ms"' 


Set 1 


3 


10 


3 


40.1% 


25.0%) 


0.0% 


Set 2 


3 


5 


3 


45.3%) 


43.0%o 


0.7% 


Set 3 


3 


3 


3 


48.3%o 


54.7%o 


7.3% 


Set 4 


4 


2 


5 


50.0%) 


64.0%) 


20.5%o 


Set r 


3 


10 


3 


39.5% 


14.8% 


0.0%. 
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three conditions we get a higher percentage of acceptable results. 
Table |4] includes eight different examples. First four use the set 
of criteria described above and for the rest we used the modi- 
fied second criterion: \a)\- (jj'^\> Bcr^'^-i^'^. The results for that 
set are denoted by an apostrophe. The case described above is 
Set 1. It is worth noticing that both criteria give very similar 
results but from the mathematical point of view (i>\ and ti;^ may 
be better indicators of ability to measure relativistic effects. 

In order to increase our chances of determining the inclina- 
tion just from an RV data set, we need to focus on eccentric bina- 
ries. With circular orbits it is not possible at all. Higher eccentric- 
ity means higher contribution of relativistic effects. With a given 
binary we may try to increase the number of RV measurements 
above 30, constant in our set, and to choose moments when spec- 
tra are taken. With the numerical model we may easily determine 
where the deviation from the Keplerian orbit is the highest, and by 
this increase the efficiency of observations. 



5 APSIDAL PRECESSION AND RV ACCURACY 

The apsidal precession is an important parameter, having great im- 
pact on the RV curve. The main cause of this effect may be tidal or 
relativistic precession. A misalignment between the star's spin and 
the orbital angular momentum of an e clipsing binary star or a third 
body may cause precession as well jMazehluOOSh . To precisely 
measure the orbital and physical parameters of binaries we have 
to take into account this effect. In the case of eclipsing binaries, it 
may also be used to confront the models of stellar structu re with 
obser vations via the stellar radial concentration parameter dMazehl 
I2OO8I) . Three examples of dw/df = 10"', \Q-\ 10"* rad yr"' are 
presented in the left-hand column of Figure |6] In our simulations 
we used units of radians per year to describe the effect. This ap- 
proach is insensitive to the different periods of systems used in the 
simulations. The pool of synthetic binaries consists of 70529 ob- 
jects, and thirty simulated RV measurements for each component 
are randomly distributed over the four year time span. Only one in- 
clination, 93°, was chosen. The ranges of the starting values of oj 
and e are presented in Table |3] 

The results of our simulations are presented in Figures |6]and 
[T] In the first figure the data used in simulations contained apsidal 
precession, but only the Keplerian model was fitted. In the second 
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Figure 4. Two binary systems are shown. The first one is presented in the panel (a) and the other one in the panels (b) - (d). They have different eccentricities 
of and 0.2. The inclination is 29° and 40°, respectively, for the first and the second system. The common parameters for both binaries are the mass of the 
first and the second component of 1 and 0.5 Mq, <jJ = 108° and peiiod P = 5 days. The most important difference comes from a non-zero e in the second 
case. Because of this the transverse Doppler effect and the gravitational redshift change with time. The Keplerian part of the RV curve is scaled and shifted 
for clarity of the figures. The semi-amplitude of the effects is quite high, however it is not unusual. From our sample we estimate that the minimal difference 
between the Keplerian RV curve and the curve with all three additional effects to be of the order of few meters per second. The average difference is equal to 
58.62 ms"' and in the extreme case of high eccentricity binaiies, can go as high as a few hundreds meters per second. 



case we fitted an additional parameter datfu/dt and compared the 
results of the best fit values with the real value of the parameter 
dcjreat /dt. In this simulation we assumed two conditions similar to 
the ones used in Section |4] The parameters used were A = 2 and 
5= 10. 

It is important to note the threshold where the apsidal mo- 
tion can no longer be absorbed by a simple Keplerian model and 
requires an additional parameter to describe observations. It is 
w„„i/dt = 10-*, 10-\ and 10"- for ctrv = 1, 10 and 100 ms"' 
respectively. The threshold was estimated with simulations and is 
based on the systems meeting the following criteria: 

• The total RMS of fitted Keplerian model is lower than three 
times the white noise applied to the RV data. 

• The error of do) is ten times lower than da>. 

A recent paper bv lKonacki et alJjaOlOt) includes five binaries 
where the RV precision reached the level of 2-10 ms"' . Four objects 
showed no apsidal motion, so with our simulations we may say that 
dojjdt = 10"^ rad yr"' is ruled out in those systems. However one 



system, HD78418, showed a small value of the apsidal motion. It 
is (1.5 ± 0.4)10"^^ rad yr"', still on the edge of detectability. Those 
results comply with the conclusions of our simulations connecting 
the precision of RV measurements with the limits on da)/dt. 

We checked the relation between the simulated time span of 
the observations and ability to detect apsidal motion. With six dif- 
ferent time spans: 1, 2, 3, 4, 10 and 100 years, we established a 
simple relation presented in Figure |8] It is an upper limit for the 
systems that will show the apsidal motion within the time span of 
the survey. The limits are based on the group of simulated objects 
fully meeting the two criteria described earlier We decided to use 
bigger groups of results rather than a few single objects as the sim- 
ulations run on a non-continuous set of input parameters. 



6 CONCLUSIONS 

Modem precision of RV measurements of binary stars reached a 
point where the subtle effects coming from the tidal distortion, rel- 
ativistic effects and apsidal precession can no longer be ignored. 
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Figure 5. Two figures showing all the systems where the inclination could be measured, based on thirty simulated RVs for each component of a binary star. The 
measurement points are randomly chosen from the whole orbital phase. The eclipses are avoided. All three conditions apply, according to the Set 1 parameters 
from Table|4] There are no results for ctrv = 100 ms"'. Two sets of binaries are distinguished, one (with inclinations lower than 20 degrees) allows much 
higher precision in inclination measurement, but such systems will be hard to observe with high RV precision due to the specti'al lines blending. The second 
set(; < 90 degrees and i > 20) degrees is more likely to return valid measurements in real observing conditions. 
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Figure 8. The relation between the observing time span and detectable 
a)„alldt for the three dift'erent RV measurement precisions. The true limit 
for a particular system may be even a bit lower, as this is an upper estima- 
tion of the limit. However we expect all the systems with a>,,,aildt observed 
with sufficient RV precision to yield the correct, non-zero rfw/rff. 



The majority of synthetic close binaries we e xamined required cor- 
rections for those effects. A recent paper bv lKonacki etai] feoioh 
shows an RV precision of up to 2 ms"'. With this level of preci- 
sion, we must model the impact of tidal distortions of components 
of binary stars on the RV curves. When this effect is ignored it 
may mimic small eccentricities in otherwise circular orbits. On the 
other hand, the relativistic effects not only need to be accounted for 
at such an RV precision, but also can be used to determine the in- 
clinations of non-eclipsing binaries based on the RV data only. In 
favourable cases, the precision of i± 0.1° should be possible. Fi- 
nally, we demonstrate that with precision RVs we can easily detect 
an apsidal motion lower than d(jjldt = 10""* rad yr"' , and in this way 
contribute to the understanding of dynamics and internal structure 
of the components of binary stars. 
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Figure 6. The left column, panels (a), (c) and (e), contains examples of simple RV curves of both components with ddyjcLt = 0. When we add doi/dt = 10"', 
10"', 10"* to those three different RV curees, we can calculate the difference between data with and without an apsidal motion. In the first case the difference 
is huge, and from the shape of O-C we may deduce an apsidal motion in the system. In the panel (c) it is still visible, but in the case of (e) it may be easily 
missed and incorporated into the total white noise. In the right column, panels (b), (d) and (f), two major groups of results may be distinguished: one group 
with a high cry,, comes from the fact that a random distribution of input pai'ameters and measurements may result in a failure of the least square method. 
However, most of the systems are in the second group around the value of the white noise ctrv for low cij„ai/dt. At some point, as the tOreai/d' increases, the 
results start to drift off the constant line in the direction of a higher RMS of the fit. The colour and the size of circles denotes log <Tfi,. 
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Figure 7. A comparison of (D^aildt and lOfuldt for 1 ms"' (left) and 10 ms"' precision (right). The figure for 100 ms"' precision is not presented here. It is 
very similar to the two results shown here, but with correct results starting around ui^calldt = 0.01. It is consistent with the conclusions from Figure|5](f). Two 
conditions similar to these in Section|4]are applied to the results. The parameters A = 2 and B = 10 are conservative but guarantee correct results. The colour 
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